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Abstract 

We consider three models (elliptic, flat and hyperbolic) of Gaussian 
random analytic functions distinguished by invariance of their zeroes 
distribution. Asymptotic normality is proven for smooth functionals 
(linear statistics) of the set of zeroes. 

Introduction and the main result 

Zeroes of random polynomials and other analytic functions were studied by 
mathematicians and physicists under various assumptions on random coeffi- 
cients. One class of models introduced not long ago by Bogomolny, Bohigas 
and Leboeuf Kostlan [IB], and Shub and Smale has a remarkably 

unique unitary invariance: 

". . . indeed it has no true freedom at all. It is (statistically) unique 
in the same sense as 'the Poisson process', or 'the thermal (black 
body) electromagnetic field' are unique. . . " 
Hannay 13, p. L755] 

Following Hannay fT2"j . we use the term 'chaotic analytic zero points' (CAZP, 
for short). We consider here three CAZP models: the elliptic CAZP, the 
flat CAZP, and the hyperbolic CAZP called by Leboeuf P3 p. 654] SU(2), 
W±, and 577(1, 1) in accordance with the symmetry group of the model. 
These models may be described analytically or geometrically. The analytical 
description is short and elementary: CAZP is the (random) set of zeroes of 
such a (random) analytic function ip, 

(0.1) m = i^C k J^ ~ 1) --- J j L ~ k + 1) z k (elliptic, L= 1,2,...), 

k=o 
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(0.2) ^) = X>VTr 



(flat, L > 0), 



k=0 

oo 



(0.3) ^) = EW L(L + 1) ' 'f + k ^ h ^ L > 0). 

fc=0 

where Co? Ci» • • ■ are independent standard complex Gaussian (random) vari- 
ables; that is, the distribution A/c(0, 1) of each (k has the density 7r _1 exp(— |z| 2 
with respect to the Lebesgue measure m on C, dm(z) = (dRez) (dlmz) . For 
the elliptic CAZP one could assume equivalently that (Co, • • • , (l) is uni- 
formly distributed on the sphere |Co| 2 + • • • + |Cl| 2 = const (which changes 
the function ip without changing its zeroes). The analytic function (jO.lj) is 
a polynomial of degree L, the function ()0.2j) with probability one is an en- 
tire function (indeed, limsup |Cfc| 1//fc = 1 a.s.), and the function (jO.Hj) with 
probability one is analytic in the unit disk. 

Why just these coefficients? Because of symmetry of the models revealed 
by a geometric description, given in Sect. 1; readers not bothered by this 
question may skip that section. 

We introduce unified notations: M. for the domain of the random func- 
tion, and G for the symmetry group. 





M. (domain) 


invariant 
metric 


G (symmetries) 


Elliptic 

Flat 
Hyperbolic 


C U {oo}, that is, § 2 

c 

D = {z G C : \z\ < 1} 


\dz\ 
1 + \z\ 2 

\dz\ 

\dz\ 
l-\z\ 2 


z^ a * + b , |a| 2 + |6| 2 = l 
—bz + a 

z i— > az + b, a = 1 

z^ aZ + b , \a\ 2 \b\ 2 = l 
bz + a 



A symmetry g G G transforms the random function ip into another ran- 
dom function z \— > ip(g(z)) whose distribution depends on g. However, for 
every g G G there exists a complex- valued function u g on A4, \u g (z)\ = 1, 
such that two random functions 



(04) z^ ^ {9{Z)) 

are identically distributed. Here 



and z I— > u g (z) 



\Mz)\\ = (nmn 



2\V2 



[l + \z\ 2 ) L / 2 (elliptic), 

exp(L|z| 2 /2) (flat), 

(1 — \z\ 2 )~ L l 2 (hyperbolic). 
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For example, a shift g(z) = z+b (for the flat case); here u g (z) = exp(iL lm(zb)). 
The symmetry ensures that the distribution of the (random) set of zeroes of 
ip is invariant under G. For details, see Sect. 1. 

In each of the three cases we have a parameter L which is assumed to 
be large. Increasing L, we increase accordingly the mean number of random 
zeroes per unit area of Ai. 

In the flat case, there is another interpretation: introducing L, we just 
make a homothety of the plane with coefficient vL. This makes some com- 
putations simpler. 

In the other two models, changing the value of the parameter L, we can 
change the properties of the process. We mention a recent discovery of Peres 
and Virag [21] that in the special case L = 1 of the hyperbolic model, the 
point process is a determinantal one [2E]. I n this case, Peres and Virag found 
a simple expression for n-point correlation functions and explicitly computed 
the distribution of random zeroes in discs in D. 

Below, the asymptotic normality is stated for smooth functionals of the 
random set of zeroes, namely, for the random variables (often called the linear 
statistics of zeroes) 

(0.5) Z L {h) = h ( z ) = I hdn ^ 

z:^ L (z)=0 JM 

where h is a smooth test function, and is the counting measure on the 
set of zeroes. In the elliptic case, n^ L (E> 2 ) = L; in the other two cases 
is infinite but locally finite. The choice of the test-functions depends on 
the model: in the elliptic case h is a real-valued C 2 -function on the Riemann 
sphere, in the flat case h is a real- valued C 2 -function with a compact support 
in the plane, and in the hyperbolic case h is a C 2 -function with a compact 
support in the unit disk. Multiple zeroes may be ignored, as well as a possible 
zero at oo in the elliptic case, since it happens with zero probability. 

Expectation of the linear statistics. For the elliptic model, invari- 
ance itself gives us the expectation of the random variable Z^{K)\ in the other 
two cases it defines the expectation up to a numerical coefficient which can be 
easily found using the Edelman-Kostlan formula Kn^ = (27r)~ 1 Alog ||^|| dm 
(see i, [23): 

EZ L (h) =L~ I hdm* 
11 Jm 



where 



(flat), 

dm*(z) = t (1 + \z\ 2 )- 2 dm(z) (elliptic) 




(hyperbolic) 
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is the invariant measure on Ai. 

Fluctuations of the linear statistics. These are computed in Sec- 
tion 2. We get 

(0.6) Var Z L (h) = j\\A*h\\l 2{m ^ + o(L~ l ) , L - oo , 

where k is a numerical constant (the same for each of the three cases), and 
A* is the invariant Laplacian on M.: 

( Ah{z) (fiat), 
A*h{z) =1(1 + \z\ 2 ) 2 Ah(z) (elliptic), 

[(1 - \z\ 2 ) 2 Ah(z) (hyperbolic). 

This shows a decay of the fluctuations of the linear statistics. 1 Note 
that in the flat case formula (jO.fij) was found in a paper by Forrester and 
Honner JU] 5 i n which they also conjectured the asymptotics in the elliptic 
case and checked it numerically. The estimate Va.r Zi(h) = 0(1), L — ► oo, 
was obtained by Shiffman and Zelditch in [2*2*] . 

Main Theorem. In each of the three cases, the random variables 

(z L {h)-L-l holm*) 

VK||A*/l||£3 (m .) V n J M J 

converge in distribution to A/"(0, 1) for L — > oo. 

The proof of Main Theorem is based on the asymptotic normality the- 
orem for non-linear functionals of Gaussian processes. Results of this type 
are known; usually, their proofs rely on the classical method of moments 
combined with the diagram technique (see Breuer and Major 7\, and the 
references therein). In Section 2 we shall prove another result in that spirit 
using a similar strategy. In Section 3, we deduce Main Theorem from this 
result and prove the asymptotics IjO.fij) . 

For real random polynomials Pn(x) = Ylk=o^ kxk wnere £fc are indepen- 
dent identically distributed (real- valued) random variables such that E£ fe = 0, 
P (f fc = 0) = 0, and E|£ fc | 2+e < oo, Maslova HHJ, [H] evaluated the variance 

^^The situation changes if we allow non-smooth test functions h. In particular, for 
the flat model, the variance of the number of random zeroes in a smooth domain A is 
asymptotic to \[L times the perimeter of A, see ^01- The difference between these two 
types of behavior reflects high frequency oscillations of random zeroes which are not taken 
into account when the test function is smooth. 



4 



of the number of real zeroes of P/v, and proved the corresponding version 
of asymptotic normality. Her methods are quite different from ours. Prob- 
ably, Main Theorem persists for more general models like zeroes of random 
holomorphic sections of high powers of Hermitian line bundles over Kahler 
manifolds extensively studied by Bleher, Shiffman and Zelditch [22JI2JIHj- 

Three toy models for the flat CAZP. It is instructive to com- 
pare the flat CAZP with simpler ('toy') models of random point processes 
in the plane, especially, random perturbations of a lattice. The first toy 
model: each point of the lattice a/^Z 2 = + If) '■ k,l G Z} is 

deleted at random, independently of others, with probability 1/2; the re- 
maining points are a random set. The corresponding smooth linear statistics 
Z^Qi) = J h{z/VZ) dnF>(z) (where is the counting measure on the ran- 
dom set) is asymptotically normal with parameters KZ^\h) = Ltt^ 1 J hdm 
and VarZ| (h) = 7r _1 J h 2 {z/y/L) dm(z) = L7r _1 ||/i||^ 2 ( m ), which is quite 
different from (jO.fjjl : the variance grows, not decays; and \\h\\ appears, not 
\\Ah\\. 

The second toy model: points of the lattice -y/7rZ 2 move independently, 
forming the random set {y/ir(k + If) + cq^i : 6 Z} where rj^i are in- 
dependent standard complex Gaussian random variables, and c G (0, oo) is 
a parameter. The corresponding Z^\h) = J h(z / 'yfD) drffl '(z) is asymp- 
totically normal with parameters ¥,Z^\h) = (1 + o(l))Ln~ l J hdm and 
VaxZf\h) = (c 2 + o(l))(27r)- 1 ||V/i|| 2 :2(m) (where V stands for gradient). 
This time, the variance does not grow, but still, it does not decay; and WVh\\ 
is not ||A/i||. 

The third toy model reaches asymptotic similarity to the flat CAZP, but 
is more complicated. Lattice points are initially aggregated into three-point 
clusters, and each cluster scatters in a special (equiangular and equidistant) 
way. Namely, we consider the random set 

{V3ir{k + li) + ce 2mm/3 r] k j : k, I G Z; m = 0, 1, 2}, 

where rjk,i and before. The corresponding Z^ (h) = J h(z / V~L) drffi (z) 

is asymptotically normal with parameters MZ^\h) = (1 + o(l))L7r _1 J hdm 
and VaiZ^(h) = (1 + o(l))L- 1 3c 4 (16vr)- 1 || A/i|| 2 2(m) . This time, we can 

mimic (IQ.fiJ) by choosing c = 2(7tk/3) 1 ' 4 . 

We come to a vague idea of two conservation laws for random zeroes: 
mass conservation and center-of-mass conservation. Mass is conserved in the 
second and third toy models; center-of-mass is conserved in the third toy 
model only. Both should be conserved in the flat CAZP, in some sense (to 
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be understood). An attempt to understand the mass conservation is made 
in our next work The center-of-mass conservation remains unexplored. 

About the three directions (e 27mn / 3 : m = 0, 1, 2): four or more directions 
may be used equally well, but two directions are not enough. Indeed, every 
quadratic form Q : C -> R satisfies ^m=iQ( e2 ™ /n ) = | A Q(0) provided 
that n > 3; however, Q(l) + Q(—l) is not proportional to AQ(0). 

Maybe, the failure of two-point clusters hints at a third conservation law. 

Acknowledgment. We thank Leonid Pastur, Leonid Polterovich and Zeev 
Rudnik for useful discussions. 

1 Geometrical description of models 

By the geometrical description of CAZP we mean something like this: 

The CAZP process on a complex 1-dimensional analytic manifold 
is the intersection of its isometric embedding into a projective 
space with a random hyperplane. 

However, we do not formalize the description; we only explain, how it works 
in the three models considered. For related more advanced theories, see the 
works of Gromov fH\ Sect. 3.3], Shub and Smale |23j . and Bleher, Shiftman 
and Zelditch 

1.1 Three homogeneous spaces 

Each model is based on a simply connected homogeneous space Ai; it may 
be thought of as a real 2-dimensional manifold or a complex 1-dimensional 
analytic manifold. 2 The former may be embedded into a Euclidean space, 
the latter — into the projective space P(C 2 ) of 1-dimensional subspaces of 
C 2 . The complex plane C is embedded into P(C 2 ) by z i— ► {(u, uz) : u 6 C}; 
in order to cover the whole P(C 2 ), one additional point is needed, oo i— > 
{(0, u) : u G C}. The symmetry group G in the real case is a subgroup of 
the group of motions of the Euclidean space. In the complex case, G is a 
subgroup of the group of projective transformations of P(C 2 ). The latter is 
covered by the group SL(2), and we may take G C SL(2); such an action 
need not be effective (since (—1) acts trivially), which is harmless. 

2 The choice of a complex structure on our real manifold does not introduce arbitrari- 
ness; there exist only two invariant complex structures, conjugate to each other, both 
leading to the same CAZP. 



6 



Elliptic model. The real manifold is the sphere 

§ 2 = {(x , x u x 2 ) : x 2 + xj + xj = 1} C R 3 , 

the symmetry group being 50(3) (orientation-preserving rotations of R 3 ). 
The complex manifold is the whole P(C 2 ), with G = SU (2); or alternatively, 
CU{oo} with transformations z i— > (az+b)/(cz+d), (* ^) G SU(2), preserv- 
ing the spherical metric |g?z|/(1 + |z| 2 ). The correspondence (well-known as 
the stereographic projection) is z = (xi + 1x2)/ (I — xo). See also [TZ| Sect. 2] 
and [20, Chap. 4]. 

Flat model. The real manifold is the plane R 2 with the group of (orientation- 
preserving) motions, that is, shifts and rotations. The complex manifold is 
the complex plane C with transformations z 1— > e t{f (z + u), preserving the Eu- 
clidean metric \dz\; or alternatively, P(C 2 ) minus a single point, with trans- 
formations preserving the point. The correspondence is just z = X\ + 1x2- 

Hyperbolic model. The real manifold is {(x , ^1,^2) : x = a/1 + x\ + x%} 
C R 3 , the upper sheet of a hyperboloid (in other words, a pseudosphere) , the 
symmetry group being the connected component of 5*0(2, 1) (rotations of 
R 3 that preserve the sheet and orientation). The complex manifold is the 
disc © = {z G C : \z\ < 1}, with transformations z 1— > (az + b)/(bz + a), 
\a\ 2 — \b\ 2 = 1, preserving the hyperbolic metric — \z\ 2 ); or alterna- 

tively, the corresponding part of P(C 2 ), with G = SU(1, 1) (isomorphic to 
SL(2, R)). The correspondence is z = (x% + ix 2 )/(l + x ). See also [T71 
Sect. 2] and [13 Chap. 5]. 

1.2 Enlarging symmetry 

The symmetry group G acts transitively on M. It acts also on discrete 
subsets of J\A, but not transitively. A G-invariant random discrete subset of 
M. is far from being unique in distribution. 

The key ingredient of the construction is an embedding 1 of Ai into a 
high- dimensional projective space P(C n ), n G {1, 2, . . . } U {00}, and use of 
the high-dimensional symmetry group U(n) of rotations of C n . Of course, 
the image l(M) is not £7(n)-invariant. However, we use [/(n)-invariance for 
determining a probability measure on hyperplanes. Choosing at random a 
hyperplane, we observe its intersection with l(A4). The case n = 00 does 
involve some technicalities; we will return to the point later. 

The embedding must be (G, L r (n))-invariant in the following sense: for 
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every g G G there exists U EU (n) such that the diagram 



M 9 -^M 



P(c n ) £¥L p(c n ) 



is commutative. Surprisingly, such an embedding is unique (up to a rotation)! 
Uniqueness of i leads to a model, (statistically) unique in the same sense as 
the Poisson process (recall the quote from Hannay in the Introduction). 
Now we switch from a heuristic to rigorous style. 

Let if be a Hilbert space over C, either finite-dimensional, or infinite- 
dimensional and separable. The projective space P{H) is, by definition, 
the set of all one- dimensional subspaces of H. Each non-zero vector x G 
H \ {0} spans such a subspace P(x) G P(H). A transformation P(U) : 
P{H) — > P(H) corresponds to every unitary operator U G U(H); namely, 
P(U)(P(x)) = P{Ux). Note that P(e^U) = P{U). Two maps i, t' : M -> 
P{H) are called unitarily equivalent, if = -P(i7) o l for some G U(H). A 
map i : — > P{H) is called holomorphic, if locally (in a neighborhood of 
any point of Ai) it is the composition of some holomorphic map M. — > i/\{0} 
and the canonical projection H \ {0} — > P{H), x i— > P(x). 

The well-known Fubini-Study metric 3 on P(H), 

dist( P(x),P(y)) = arccos , , 

||x||||y|| 

is f/(n) -invariant. Given a one-to-one map t : — > P(H), we get a metric 
on A4, 

dist t (;z, z') = distu(2;), ; 

assuming smoothness of l and taking z' infinitesimally close to z we get a 
tensor field g L on A4, the Riemannian metric induced by i. If t and t' are 
unitarily equivalent, then they induce the same Riemannian metric (since 
dist t = dist t '). 

Calabi's rigidity theorem. If two holomorphic embeddings 4 of a complex 
manifold into P(H) induce the same Riemannian metric on M., then they 
are unitarily equivalent. (0 Th. 9], see also [23 Sect. 4].) 

3 Or rather, the geodesic metric corresponding to the Fubini-Study metric tensor, up to 
a coefficient. 

4 'Embedding' means one-to-one. The theorem holds also for immersions, but we need 
only embeddings. 
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If a holomorphic embedding i : M. —>■ P{H) is (G, ?7(n))-invariant (recall 
(fl.lj) ). then it induces a G-invariant Riemannian metric on M.. However, 
such a metric is unique up to a coefficient (since G can rotate M. around any 
given point). It means that the coefficient is the only possible parameter of 
a holomorphic {G, U (n)) -invariant embedding i : M. — > P(H) (treated up to 
unitary equivalence). 

Another implication of Calabi's theorem: if a holomorphic embedding 
l : M. — > P{H) induces a G-invariant Riemannian metric on M., then i is 
(G, U (n)) -invariant. 

By the way, (G, £7(n))-invariance leads to a projective representation of G 
(that is, a homomorphism from G to the factor group of U(H) by its center), 
not necessarily a homomorphism G — ► U(H). 

1.3 Relation to the analytical description 

In order to finish the description, it remains to write down the corresponding 
embeddings t in each of the three cases, to check invariance of induced met- 
rics, and to explain the rotation invariant choice of the random hyperplane 
in the case when the dimension of H is infinite. 

All invariant metrics on M. are proportional to the one mentoned in 
the Introduction. The only freedom left to us is to choose the numerical 
coefficient y/~L of the invariant metric on M. . 

Elliptic case. For the sphere § 2 , the existence of such an embedding 
depends on the parameter L. Only for L = 1,2,3,... the invariant metric 
on § 2 is embeddable into P(H). 5 For such L, an embedding S 2 — > P(C i+1 ) 
is well-known to physicists as the system of spin-J coherent states (J = 
L/2), see I2H1 Sect. 4.3], j^Zj- Sometimes, mathematicians call the wave 
function of a coherent state "the Szego kernel" (see [3]). Treated up to 
rotations of C L+1 , the map is unique. By rigidity, every such embedding uses 
only a finite-dimensional subspace of H. 

The embedding S 2 — > P(C L+1 ) may be described as the composition 
§ 2 -> C -> C L+1 \ {0} -> P(C L+1 ), where § 2 -> C is the stereographic 
projection mentioned in Sect. 1.1, C L+1 \ {0} — ► P(C L+1 ) is the canonical 
projection x i— > P(x), and t : C — >• C L+1 \ {0} is given by 




which evidently corresponds to (|0.1|) . 

5 One can easily get this using the Edelman-Kostlan formula: if L is non-integer, then 
we would get a non-integer answer for the average number of zeroes of our random function 
on § 2 . 
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The induced metric is easy to calculate: 

L 

(Z(z), l{z')) = E V© ^ V© ^ = (1 + rtf 

k=0 

distt^, 2 J = arccos — — = arccos 



v/|2 



1^)1111^)11 VxATRyiTk 

G-invariance of the metric for every L follows immediately from its G-invariance 
for L — 1; however, for L = 1 the embedding C — > P(C L+1 ), given by 
z i— > : w G C}, is already familiar to us (recall the beginning of 

Sect. 1.1). An explicit calculation gives for Az — > 

V ( 1 A 1 A 

dist t ( 2 — ^Az, z + 



arccos I 1 — 





IA2 


2 


(1 + 1 




4l 


A*| : 


2 ) 2 - Re 2 (zA^) 



L/2 



( L \Az\ 2 , . . \ / - 
arccos ^1 - - (1 + | ^ |2)2 (1 + o(l)) 1 = VL- 









1 + 





2 



(1 + 0(1)). 



Flat case. If the needed embedding 1 : C — ► P(iP) exists for L = 1, 
then for every L G (0, 00) the map z 1— > i{L x l 2 z) fits; we restrict ourselves 
to L = 1. The construction of t is well-known to physicists as the usual 
system of coherent states (Schrodinger, von Neumann, Klauder, et al.), see 
|2lH Chapter 1]. Its explicit form is given (similarly to the elliptic case) by 
the composition C — > H \ {0} — > P(H) where the first map Z : C — > H \ {0} 
is given by 



z 2 z 3 



Z(z)= [l,z,— I el = H 



2!' a/3! 



l2 



cf. (jnn. 

We have 



00 y. 

iZ(z), Z(z')) = V -7= -7= = exp(z^) 
7—: v fe! v fc! 

fc=0 



dist t (z, z') = arccos exp ^Re^z') — — | | 2 — - 



1 2? 



X i /12 
arccos exp ^ — - 1 z — z \ 

dist L (z,z + Az) = arccos ( 1 - -|A^| 2 (1 + o(l)) J = |A*|(1 + o(l)) . 



2' 
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Hyperbolic case. The corresponding system of coherent states is usually 
investigated in terms of unitary representations of the group G = SU(1, 1), 
see J201 Chap. 5], [13 Sect. 2], which leads to a special treatment of integer 
values of L. In our approach, projective representations appear irrespective 
of unitary representations, and L runs over (0, oo). 

As before, l(z) = P(Z(z)), but now z G D, and I : D — > I 2 is given by 



M = (l, VLz, 4^z\ y/MtdW&f, . . .) , 

cf. ()U.3|) . The corresponding embedding of the hyperbolic plane to P{H) was 
known already to Bieberbach in 1932. We have 



\1 — zz: 



z' k = (l-zz')- L 



fc=0 

dist t (2, z') = arccos 



IzlVl - Iz'l 2 



arccos 





|A2 


2 


(1-1 


z\ 


2 1| 
4l 




2 ) 2 -Re 2 (^Az) 



-L/2 



arccos (l - | JM^ (1 + o(l))^) = ^^^(l + o(l)) . 













2; 


2 



1.4 Random hyperplane in Hilbert space? 

In C n , the only [/(n) -invariant distribution on hyperplanes, the uniform dis- 
tribution, may be represented via the normal (to the hyperplane) vector 
(£1, . . . , Cn) distributed uniformly on the sphere. The normal distribution for 
(Ci, • • • , Cn) can be used instead (as well as any spherically invariant distri- 
bution). 

In a Hilbert space H = I 2 there is no rotation-invariant probability mea- 
sure on hyperplanes (nor Lebesgue measure in H). Nevertheless, an infinite 
sequence of independent normal Ac(0, l) random variables Cfc can be used. 
Of course, it is not a random element of H, since the event I Cfc 1 2 < 00 is 
of probability 0. However, for every x = (ci, c 2 , . . . ) G I 2 the series ^CkCk 
converges a.s., and we may denote its sum by (x,Q. The 'bad' set of zero 
probability, on which the series does not converge, depends on x; the union 
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over all x G H is not a set of zero probability. 6 We cannot choose £ at 
random and speak about 'the function x i— ► (x, on if'. 

What we really need, is the function z \— > (£(2), C) f° r 2 £ that is, the 
function x 1— > (x, for x G ?(-M) only. The corresponding series converges 
a.s. for all these x simultaneously. Thus, 'the random hyperplane of W is 
ill-defined, but still, its intersection with 1{M) is well-defined. 

It remains to explain unitary invariance of the construction described 
above. Random variables (x, £) can be used simultaneously for a countable 
set of points x. In particular, for every orthonormal basis (ei,e2, ...) of 
H, the sequence of random variables (e^, C) is well-defined. Some reflection 
shows that the joint distribution of these (e&, C) does not differ from that of 

It follows easily that our construction can start with an arbitrary basis, 
resulting in the same distribution of the random function z *—>■ (l(z), (). See 
also [14, Sect. 1.3]. 

2 Asymptotic normality for non-linear func- 
tional of Gaussian processes 

2.1 The result 

Let T be a measure space endowed with a finite positive measure /1, fiiT) = 1. 
A complex-valued Gaussian process on T may be defined as 

(2.1) w(u,t) = J2Ck(u)fk(t) 

k 

where fk : T — > C are measurable functions such that 

^|/ fc (t)| 2 <oo foralliGT, 
k 

and (k — 6c + if]k are independent standard complex Gaussian variables; i.e. 
Cfc ~ Ac(0, 1). The latter means that and % are independent centered 
Gaussian (real) variables with variance |; then = and E|(^| 2 = 1. We 
restrict ourselves to the case 

Y^\fk{t)\ 2 = l for all* GT. 

Then tu(*) ~ W c (0, 1) for all t G T. 

6 In fact, it is the event ^ \(k\ 2 = 00, of probability 1. 
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Now, a few words about the convergence of the series (j2.1j) . Treating 
each term ( k (u)f k (t) as an element of the space C 2 = L 2 ((f2,P) x (T,p)), 
we have 

HCfc/fclU 2 — IICfc||L 2 (0,P)||/fc||L2(T, At ) = ||/fc||L 2 (7» 

and the terms are pairwise orthogonal. Therefore, the series (j2.1j) converges 
in the space 7 C 2 . 

In what follows, we always assume that the sum of the series (|2.1|) is not 
the zero function in the space C 2 . 

The correlation function p : T x T — > C of the process w (t) equals 

P (s, t) ^ e {w(s)W)} = e f^y^) ■ 

k 

Clearly, \p(s,t) \ < 1 and p(t,t) = 1. 

Consider a sequence of complex Gaussian processes w n with the correla- 
tion functions p n (s,t) and define a sequence of random variables 



Z n = J <p(\w n (t)\)Q(t)dn(t) 
where (p : M + — > R is a measurable function such that 

POO 

/ (p 2 (r)e~ r l 2 rdr < oo, 
Jo 

and G : T — ► M. is a measurable bounded function. We shall prove that 
under some natural assumptions on the off-diagonal decay of the correlation 
functions p n (s,t) when n — > oo, the random variables Z n are asymptotically 
normal. 

2.2 Theorem. Suppose that for each a G N 

, 0Q , . ( JJ T2 \Pn(s,t)\ 2a Q(s)Q(t)dp(s)dp(t) 

(2.3) hmmf — i p- — > 0, 

n-oo SUp sgT J T |p n (s,t)|d/i(t) 

and i/iai 

(2.4) lim sup / \p n (s,t)\ dp{t) = 0. 



7 Moreover, for each t € T, the series H2.lt converges in L 2 (Q, P) (for an obvious reason), 
and by the Kolmogorov-Khinchin theorem ^3 Ch. 3] the series converges in L 2 (T,/j,) for 
almost all to E O. 
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Then the distributions of the random variables 



Z n — ~E,Z r , 



VVar Z n 

converge weakly to 7V(0, 1) for n — > oo. 

// (p is an increasing function, then it suffices to assume that condition 
()2.3j) holds only for a — 1. 



Remarks: 

(i) The function \p(s,t)\ 2a which appears in the integrand on the left-hand 
side of condition ()2.3j) is a positive definite function on TxT (see formula (|2.fij) 
below). Hence, for a G N, the integral ff T2 \p(s,t)\ 2a Q(s)Q(t) dp(s)dp(t) is 
always non-negative. 

(ii) The role of condition (|2.4jl is to guarantee that 

lim Var Z n = 0. 

In fact, under assumption (|2.3jl . the sequences 

n i— > VarZ„, and n ^ J \Pn( s ^)\ dp(t) 
have the same decay (see expression (|2.7jl below). 



2.2 Some preliminaries 



We shall deal with the space = f L 2 (C, Ac(0, 1)) and its subspaces V m 
which consist of the polynomials in ( and ( of degrees a and (3 correspond- 

"2 
J C 



ingly, a + f3 = m. The space L£ has a polynomial basis 



where : £ a C : is the orthogonal projection of the polynomial ('"C onto the 
subspace if :m: = f V m Q V m -u m = a + j3 [HJ Example 3. 32]. 8 Thus any 
square integrable function $ of a random variable ( ~ Ac(0, 1) is of the form 



*(C)= E -mm ■■<??■■• »*II 2 = E 



-a/3 1 



3 The expression : C a C^ : is called the Wick product. 
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In what follows, we shall deal only with the radial functions <£>(|C|)- I n 
this case, 



oo 

C 2a IAl2a 



(2-5) ^(KD = E^ : ICI 



0=0 



where C2 a E K and J2 a C L = IMP- Indeed, the group of rotations ( 1— > e i6 *C 
acts on the space leaving invariant the subspaces H' m: . If m is odd, then 
i7 :m; contains no radial functions; if m = 2a is even, then the subspace of 
radial polynomials in H :m: is one-dimensional and is spanned by : \(\ 2a : = 
: :• 

In the proof of Theorem I2.2| we may assume without loss of generality 
that 

%(|C|)e- K|2/2 rfm(C) = 0, 



that is, EZ n = 0. Then expansion ()2.5j) starts with a = > 1, C2 ao 7^ 0. 
Since 

(2.6) E (— : \ W {s)r \w(t)\* ^ - 1 '^'^ * * = & 



a! / \/3\ / I otherwise 

[TU Theorem 3.9] (see also Example at the end of the next subsection), we 
have 

EM\ Wn (s)\M\ Wn (t)\)}= 4a\Pn(s,t)\ 2a , 

and 

/ r \ 2 

72 



(2.7) EZ^ = 2( / ^(|«; B (t)|)e(t)d/*(t) 

E{<p(\w n (s)\M\w n (t)\)}Q(s)e(t) d/i(s)dii(t) 

= J2 C 1* [[ \Pn(s,t)\ 2a Q(s)&(t)d^s)d^t) . 

The proof of Theorem 12 .21 uses the classical method of moments, though it 
cannot be applied directly to the random variables Z n since ip (and therefore, 
Z n ) need not have more than two moments. First, we approximate <p by the 
polynomials 

^(ici) = E^7 : Ki 2Q: 

z — ' a! 
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and Z n by the random variables 

Z n , m = i ipm{\w n {t)\)Q{t) dfj,(t) . 



T 



Then the moment method will be applied to Z njTn . 

Applying formula ()2.7|) to Z n — Z n , m (that is, replacing the function ip by 

(p - (p m ), we get 

E(Z n -Z n , m ) 2 = J2 °laff \Pn(s,t)\ 2a e(s)&(t)dfi( S )dfi(t) 



a>m+l 



< llftll 2 

— II w llL°°(T,/i) 



T 2 



V 4, sup / \p n (s,t)\dfi(t) 



a>m+l 

^ ll<5>ll W,M) „,„ ,,2 / / , n |2,.„. 



- J a A ll^-^llig I IJPn(s,t)r°e(s)e(t)dp(s)dp(t) 



T 2 



where 77(a) denotes the expression on the left-hand side of (|2.Hj) . Using again 
formula ()2.7|) . we finally get 

IIBII 2 

TP I <V V \ 2 ^ H w H£°°Cr» || || 2 ^,^2 

JH(Z n -Z nim J < — 2 7 — r- _ Emilia lkZ n . 

Hence 

\ 2 



r. 



lim sup E — ^=^= ' ] =0 

m ^°° n y v"VarZ n ^/VarZ n)r 



and if for each fixed m > ao the random variables Z n>m / y/Var Z nim con- 
verge in distribution to A/"(0, 1) when n — ► 00, then the random variables 
Z~ n / \/Var Z n have the same property. 

From now on, we always assume that <p is a polynomial (and hence Z n 
has the moments of all orders). 

2.3 More preliminaries (the diagram formula) 



The next step is to evaluate the moments Ez" p , p e N, and to compare them 

h (Ezir 

We have 



with {E,ZTf^ 2 times the moments of the standard normal distribution 



(2.8) ezi=!...( E|n^K(t,)i)lne(t,)^) 

9 e ^^/■••I E {n=w^ ! }n e (^ 

a 1) ...,a p >a 1 P J ^ J U=l J j=l 
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We compute the integrand E|f|^ =1 : \w n (tj)\ 2aj : j using the diagram tech- 
nique [HI Chapter 3]. 

Fix the exponents ai, . . . , a p . A diagram 7 is a graph with 2(ai + - • ■ + a p ) 
vertices labeled by the indices 1, 1, 2, 2, . . . ,p, p (atj vertices are labeled by j 
and other otj vertices are labeled by j), and each vertex has degree one (i.e., 
the edges have no common end points). The edges may connect only the 
vertices labeled by i and j with i ^ j. The set of all such graphs is denoted 
by r(ai, . . . , a p ). For some choices of ax, ... , a p , T(a±, . . . , ot p ) may be the 
empty set, for example, T (0.1,0.2) — iff o± ^ 02, and T(ai,a 2 ,as) = if 
Qi > a 2 + o 3 . The value of the diagram 7 equals 

Vy{ti, . . . ,t p ) = Yl P(U,tj) 

(i,j)£e(y) 

where the product is taken over all edges of 7. In this notation, the diagram 
formula [Til Theorem 3.12] reads 9 

(2.9) E|n:Wt,)|^:U £ V^h, . . . ,t p ) 

li=i ) 7er(ai,...,ap) 

(as usual, summation over the empty set means that the sum has zero value). 

2.10 Example. Consider the simplest case p — 2. If ai ^ then 
r(«i,a 2 ) = 0, and E{: \w(h)\ 2ai :: \w(t 2 )\ 2a2 :} = 0. Now, suppose ai = 
o 2 = o. Then we can glue together the vertices labeled by the same indices. 
We get a graph with four vertices and two edges of multiplicity a: 



(2.11) 




The edges connecting 1 and 2 contribute the factor p(ti,t 2 ) a , the other a 
edges contribute the factor p(t 2 ,£i) a = p(ti,h) a - Thus the value of the dia- 
gram is \ P (hM)\ 2a , and E{: Hti)| 2a :: \w(t 2 )\ 2a :} = $T(o,o) ■ \p(h,t 2 )\ 2 <*. 
It remains to find the total number of diagrams in T(a, a). 

All the diagrams in T(a, a) can be obtained from the fixed one by permu- 
tation of a vertices labeled by 1, and by another independent permutation of 
a vertices labeled by 2. Therefore, §T(a, a) = (a!) 2 , and we recover formula 

9 The diagrams we consider do not contain edges which join the vertices labeled by i 
and j (and by i and j) since E {w(ti)w(tj)} =¥,< w(U) w(tj) > = 0. 
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2.4 The main argument 

For p even, there are diagrams with a simple structure whose total contri- 
bution to EZP equals {EZl) p/2 E£ p , where £ ~ W(0, 1). Contribution of the 
other diagrams is negligible. 

2.12 Definition. A diagram 7 is called regular if the set {1, 2, ... ,p} is split 
into q = p/2 pairs and there are no edges between the vertices i and j if i 
and j belong to different pairs. Otherwise, the diagram is called irregular. 

In other words, the diagram is regular if after glueing together the vertices 
labeled by the same index it becomes a disjoint union of q = p/2 "elementary 
diagrams" drawn on ()2.11|) (having, generally speaking, different multiplici- 
ties of the edges). 




regular diagram 




irregular diagram 



Suppose p is even and the diagram 7 is regular. For each j, we glue 
together all the vertices labeled by j and j, obtaining a reduced diagram 7*; 
i.e., a graph with p vertices and multiple edges which is split into q = p/2 
connected components. 

For example, if 7 is related to the partition 

(2.13) {1, 2, . . . ,p} = {1, 2} U {3, 4} U • • • U {p - l,p} , 

then the reduced diagram 7* is 

2/3i 2/3 2 2/3 q 



where 2(3k = 2a 2 k-\ = 2a2k is the multiplicity of the edge which connects 
the vertices 2k — 1 and 2k. The k-th component of the reduced diagram 
7* contributes by the factor \p(t 2 k-i, ^2fc)| 2/3fe , so that the value of diagram 
7 is Ilt=i\p(hk~i,h k )\ 2 ^. 
llf=i ®(ti)dn(ti), we get 



Integrating this over T p against the measure 



(2.14) 



n 

k=l 



T 2 



\p{s,t)\ 2 ^Q{s)Q{t) dfj,{s)dfj,{t) . 



Now we need to do some counting. Each reduced diagram 7* is defined by 
the choice of the partition, like (|2.13)1 . and the multiplicities f3\, . . . , (3 q of the 
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edges. Since the glueing procedure is not one-to-one, each reduced diagram 
has a "multiplicity" (that is, the number of regular diagrams which give the 
same reduced diagram 7* after glueing) depending on fli , . . . , /3 q which also 
should be taken into account. 

Notice that if we started with a different partition of the set {1, 2, . . . ,p} 
into disjoint pairs, then anyway we would finish with the same answer (j2.14|) . 
Therefore, we must multiply the expression ()2.14|) by the number of partitions 
of the set {1,2, ... ,p} into disjoint pairs. The number of such partitions 
equals E£ p where f ~ W(0, 1) [11 Remark 1.29]. 

Next, let us count the "multiplicity" of the reduced diagram. If we fix a 
regular diagram 7, then all other regular diagrams having the same reduced 
diagram 7* can be obtained from 7 by p independent permutations: we 
can permute a± vertices labeled by 1, then a 2 vertices labeled by 2 and so 
on. Since these permutations are independent, the total multiplicity of the 
reduced diagram equals 

(2.15) aM- ■.■a p \ = ({3 1 \) 2 (f3 2 \) 2 ...((3 q \) 2 , 

and we need to put this factor before the product ()2.14|) when we sum- 
mate over reduced diagrams; i.e. over all possible choices of the numbers 

(3l,...,(3 q > Oq. 

Combining all pieces, we get 

2 2 

(regular) = £ '^'^2 • • • (A,!) 2 (product (EH) 

ft,...,0,>ao ' ' ' W 

= (W) ( £ c\J\ \p(s,t)\^Q(s)Q(t)d^s)d^t)) W (EO (EZl) p/2 . 
Therefore, 

W n = (EO {EZ 2 n ) p/2 + (irregular) , 

and we need to show that the contribution of irregular diagrams is negligible 
with respect to the main term (1EZ 2 ) P ^ 2 when n — ► 00. 

Since tp is a polynomial, there are only finitely many irregular diagrams 
which enter expression ()2.8j) for EZ^ after plugging in the diagram formula 
(|2.9jl . We shall show that if the diagram 7 is irregular, then 
(2-16) ^ 

/••■/ \Vy(ti,t2,... : t p )\dfj l (tx)...dfj l (tp) = I (sup / \p n (s,t)\dn(t)) 

J J TP \ VseT Jt J 
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for n 



and 



oo. Then due to assumption (|2.3jl (and expression ()2.7|) for EZ%) 
V^tu t 2 ,..., tp) | dp{t x ) . . . dfi(t p ) = o ((E^) p/2 ' 



TP 



(irregular) = o ( (EZ^) 



n 



oo 



which finishes the proof of Theorem 1.2 in the case of general radial functions 
up to the proof of (|2.1fij) . 

Proof of (|2.1(jj) : First, we make another reduction of the diagram and define a 
p- vertex graph with simple edges which couple the vertices % and j if and only 
if at least one of the pairs (i, j) or (J, i) was coupled in the original diagrams 
(without taking into account the multiplicities of the original coupling). We 
denote the reduced diagram by 7**. For example: 





Then 



(2.17) 



a diagram, 7 



|V 7 (t! 



3 4 

its reduction, 7* 



< n w 

(i,i)ee(7**) 



1, tj, 



where the product is taken over all edges of 7**. We have to estimate from 
above the integral of \Vy\ over T p . Replacing \ Vy\ by its upper bound f)2 . 1 71) . 
we obtain the integral which factorizes into the product of integrals described 
by connected components of the diagram 7**. 

Let us start with one m- vertex component of the diagram 7**. The com- 
ponent can be a complicated graph - - anyway, we can always turn this 
graph into a tree with m vertices by deleting some edges (this procedure 
only increases the integral we are estimating). Having a tree, we choose a 
vertex belonging to only one edge and integrate it out, which gives the factor 
(sup sgT J T \p n (s,t) \ dfj,(t)) and the rest of the tree which is a new tree with 
m — 1 vertices. By induction, any m-vertex tree describes the integral which 
does not exceed 



sup 



\p n (s,t)\d(x(t) 



T 



m—l 



Now, suppose the reduced diagram 7** has k connected components and 
the i-th component has rrii vertices. 10 Then the right-hand side of ()2.17j) 



10 



Observe that mi + • ■ ■ + = p. 
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integrated over T p does not exceed 

r \ (mi-l)H h(m fc -l) / \ p-fc 

sup / \ Pn (s,t)\d/i(t) = (sup / |p n (s,t)| I 

seT Jt J \seT Jt J 

Since the diagram 7 is irregular, k < p/2 and we get (j2.16|) . 
2.5 The last step 

It remains to explain why in the case when the function (p is increasing, it 
suffices to assume that condition (|2.3J1 holds only for a — 1. 

In the proof given above, condition ()2.3j) was used only in the estimate 



SUp / \p n (s,t)\ dfi(t) < 

sgt Jt 



where a > 1 is the minimal positive integer such that c 2ao 7^ 0, and i](a) 
is the left-hand side of (12. 3 J) : i.e., we need condition ()2.3|1 only with a = oq. 
If the function (p increases, then cto = 1. Indeed, it is easy to find that 
: ICI 2 := |C| 2 -1- Then 

poo roc 

c 2 = <p(r)(r 2 - l)re~ r2/2 dr = (ip(r) - <p(l)) (r 2 - l) re~ r2/2 dr > , 



'0 jo 
concluding the proof of Theorem 1.2. 



3 Asymptotic normality for chaotic analytic 
zero points 

3.1 From chaotic analytic zeroes to non-linear func- 
tionals of Gaussian processes 

Recall the random objects defined in the Introduction: the Gaussian ana- 
lytic functions i/)l{z) ( see (EH), (|0-2j) and (|0.3jl ). the random measures n^ L 
counting their zeroes, and the linear statistics Z^ih). Since 

2itdn^ = A log dm = A* log \ip\ dm* 

where dm* and A* are the invariant measure and the invariant Laplacian on 
A4, we have 

Z L (h) = ±- [ \og\Mz)\(^h)(z)dm*(z) 
^ Jm 
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obtaining in each of the three family of non-linear functional of 

the complex- valued Gaussian process ipL defined on M.. We normalize the 
process ip putting 

[Mz){l + \z\ 2 )- L/2 (elliptic), 
(3-1) w L (z) = - 7 ^= = Mz)e-W (flat), 

V^IWJI y L ( z )(l - \ z \^)L/2 ( hy p er bolic). 

In the flat and hyperbolic cases we use a large real parameter instead of the 
large integer parameter used in Section 2. Note that the (non-Gaussian) 
random process log |w(z)| is stationary (that is, its distribution is invariant 
with respect to the corresponding group of isometries), though the Gaussian 
process w(z) is not (because of the phase multipliers in (J0.4j0 . In all three 
cases, w(z) ~ Ac(0, 1) for each z. We obtain 

Z L (h)-EZ L (h) = ^- [ \og\w L \A*hdm* . 
2tt Jm 

Indeed, the left-hand side has zero expectations since Elog = const, 

and the invariant Laplacians of the test functions are orthogonal to the con- 
stants. It remains to check conditions (|2.3J) (with a = 1) and (|2.4J) of The- 
orem I2.2| and to compute the asymptotics of the variance of the random 
variables Zl(H). 



3.2 Checking conditions (O) and RIM 

Let 



/ \ E (iIjl(zi)iPl(z2, 

Pl(zi, z 2 ) = E [wl(zi)w l (z 2 )J 



y/E\M*l)\ 2 -nMz2)\ 2 ' 

The function of two variables, \pl(zi, z 2 )\, reduces to a function of one vari- 
able, \p L (z,0)\, by G-invariance, \pl(zi,z 2 )\ = \pl{q{zi), g(z 2 ))\. Explicit 
formulas for p L (z,0) = \p L (z,0)\ follow from flQJ, JQ2J|, and (ITU : 

namely, 

p L (z,0) = (1 + \z\ 2 )~ L/2 (elliptic), 
Pl{z,0) = exp(~L\z\ 2 ) (flat), 
p L (z,0) = (1 - \z\ 2 ) L/2 (hyperbolic). 

A straightforward inspection shows that in all three cases the functions 
j^pL^z, 0) converge weakly to the point mass at the origin: 



hm I p L (z,0)e(z)dm*(z) = 9(0) 

' M 



2ir 
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for any continuous test-function (as usual, with a compact support if Ai 
is non-compact). 

By G-invariance, for any (3 > 0, any continuous test-function and any 
z 2 e M, 

lim ^ / \ PL (z 1 ,z 2 )fQ(z 1 )dm*(z 1 ) = Q(z 2 ) 

(recall that \pl\^ = \P/3l\ )• At last, multiplying the both sides of this equa- 
tion by 0(^2) and integrating by z 2 , we get 

(3.2) lim ^- ff |pi(«i,^)|"e(«i)e(ft)dm*(zi)dm*(^) = ||0||i 2(w *). 

L^oo 2n J J M 2 K 

Now, conditions ()2.3|) and ()2.4|) become evident: for L — > 00, 

2tt 



SU P / \pL(z 1 ,z 2 )\dm*(z 1 ) = / Ipl^, 0| dm*{z 



z 2 €MJM Jm 



L ' 



this gives us ()2.4|) . To get ()2.3|) . observe that the double integral in the 
numerator of ()2.3j) has the same decay 



7T 

|PL(^l,^2)| 2 0(^l)0(^2)^*(^l)^*(22) ~ J ■ ||©||za(t»*)- 

M 2 



3.3 Asymptotics of the variance 

We shall use formula (|2.7|) for the variance of Z(h), 



EZ 2 L {h) = 7^5^/7 |PL(^i^2)r cl A*/ i (^ 1 )A*/i(^)rfm*(^ 1 )^(^) 



47r2 "I ■/ ^ 2 
where C2 a are defined by the expansion 



log|C| = X^^f : ICI 2 " : forC~A/- c (0,l) 



a=0 



(cf. (Q). Denoting 

a>l 

we reduce fl0.6j) to relation (|3.2j) with /3 = 2a: 
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